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We present a theoretical description of the radio-frequency (rf) spectroscopy of fermionic atomic gases, 
based on the second-order response theory at finite temperature. This approach takes into account the 
energy resolution due to the envelope of the rf pulse. For a noninteracting final state, the momentum- and 
energy-resolved rf intensity depends on the fermion spectral function and pulse envelope. The contributions 
due to interactions in the final state can be classified by means of diagrams. Using this formalism, as well as 
the local density approximation in two and three dimensions, we study the interplay of inhomogeneities 
and Hartree energy in forming the line shape of the rf signal. We show that the effects of inhomogeneities 
can be minimized by taking advantage of interactions in the final state, and we discuss the most relevant 
final-state effects at low temperature and density, in particular the effect of a finite lifetime. 
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I. INTRODUCTION 

In many-fermion systems, the low-energy properties are 
often determined by single-particle excitations across the 
Fermi surface. The character of these excitations depends 
on the nature of the ground state, which itself depends on 
the interactions. The study of single-particle excitations 
is therefore a key to understanding the ground state and 
the role of interactions. In superconductors, for instance, 
a gap in the single-particle excitation spectrum reveals the 
condensation of Cooper pairs in the ground state. In a large 
class of materials, the interactions bring only quantitative 
changes with respect to a noninteracting ground state. The 
single-particle excitations are then similar to uncorrelated 
particles, albeit with a renormalized mass and a finite life¬ 
time. The collection of these “quasiparticles” forms a Fermi 
liquid, which can be characterized by a small number of ef¬ 
fective parameters.^ For strong interactions and/or reduced 
dimensionality, qualitative changes may occur in the ground 
state, leading to the disappearance of the quasiparticles 
and the emergence of more complex, sometimes mysterious, 
excitations.^’^ The absence of quasiparticles in a fermion 
system is a hallmark of non-Fermi liquid physics, indicating 
an unconventional ground state. 

For electronic materials, angle-resolved photoemission 
spectroscopy (ARPES) gives access to the single-particle 
excitations and allows one to probe the existence of 
quasiparticles.The signature of quasiparticles is a peak at 
low energy in the spectral function, which is the momentum- 
energy distribution of the single-particle excitations, de¬ 
noted A(k, e). Conversely, a structureless spectral function 
signals the absence of quasiparticles. ARPES experiments 
require clean surfaces and ultrahigh vacuum, and an en¬ 
ergy resolution below the typical excitation energy of the 
quasiparticles. Steady improvements in recent years and 
the development of laser ARPES have made it possible to 
measure the spectral function with excellent resolution in 
several condensed-matter systems.When it is present, the 
quasiparticle peak and its dispersion anomalies can help in 
identifying the interactions that determine the quasiparticle 


dynamics. 

Fermionic cold-atom gases open new avenues in the study 
of quasiparticles, especially thanks to the possibility of tun¬ 
ing both the dimensionality and the strength of interac¬ 
tions. Radio-frequency (rf) spectroscopy is presently the 
best method to measure the spectral function of cold-atom 
systems. Unlike in conventional ARPES, photoemission spec¬ 
troscopy in ultracold atoms is performed using rf photons, 
which carry negligible momentum but only supply an energy 
hv. The momentum of the extracted atoms is then measured 
using the time-of-flight technique. If the particles are decou¬ 
pled in the final state, their energy and momentum distri¬ 
butions follow the spectral function of the photon-induced 
hole, which is the occupied part of the spectral function, i.e., 
A(k, £)f(£), where /(e) is the Fermi function.®’® 

The interpretation of photoemission and rf experiments 
may be complicated by the unavoidable interaction in the 
final state, as well as several other difficulties. In ARPES, 
these are, for instance, the sample surface, which breaks in¬ 
version symmetry and produces interference, or the screen¬ 
ing of the electromagnetic field, which prevents light from 
entering the bulk of the material. In rf spectroscopy of 
cold atoms, the main concern is the inhomogeneity of har¬ 
monically trapped gases. When interactions and excitation 
energies are not too low, as in the studies of the BCS-BEC 
crossover,^® some of these difficulties may turn out to be 
irrelevant. For weak interactions, however, they will eventu¬ 
ally become important. If the signal is broadened by final- 
state effects, averaging over inhomogeneities, and finite 
energy resolution, a precise modeling is necessary in order 
to recover the crucial information about the quasiparticles. 

In the established theory of rf spectroscopy, one com¬ 
putes the instantaneous transition rate to the final state. 
This can be done either by linear response,® which pro¬ 
vides the current Nf of particles transferred to the final 
state, or by time-dependent perturbation theory (Fermi 
golden rule).® At leading order, Njr is related to a response 
function, which can be represented by bubblelike Feyn¬ 
man diagrams.In this approach, the effect of inhomo¬ 
geneities has been investigated at the mean-field level^^ 
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or using the local-density approximation (LDA).®’^^ To cir¬ 
cumvent the difficulties raised by inhomogeneity, a Raman 
local spectroscopy was proposed theoretically,® while a 
tomographic technique and a method to selectively ad¬ 
dress the cloud center^^ were demonstrated. Final-state 
effects have been treated in the mean-field approximation^® 
by sum-rule arguments,^^’^® within a reduced basis,a 
1/iV expansion,^® diagrammatically,^^’^^ or through self- 
consistency requirements.^^ Most of these studies have fo¬ 
cussed on the BCS-BEC crossover problem. 

For intermediate or weak interactions, the finite energy 
resolution must be considered. The relevant quantity to 
calculate is no longer Nf, but the total population Nf of 
the final state, created over the duration of the rf pulse. 
Momentum-resolved rf experiments indeed measure the mo¬ 
mentum distribution n^(fc, t) at a time t after the extinction 
of the rf pulse. If atoms were excited at a constant rate, Nf 
and Nj: would carry the same information, but this is not the 
case in practice. In this paper, we present the calculation of 
n^(fc, t) within equilibrium response theory. The derivation 
is performed in the finite-temperature Matsubara frame¬ 
work. Unlike Nf,nf vanishes at first order in the atom-light 
coupling. At second order, the momentum distribution is 
related to a three-point response function, whose contri¬ 
butions can be classified using Feynman diagrams. These 
diagrams have three external vertices, unlike the bubble 
diagrams of the established theory, which have only two. 
The leading contribution reproduces the known result,®’® 
albeit convolved with a resolution function, which depends 
on the envelope of the rf pulse and on the spectral function 
in the final state. 

Simulations based on this formalism have been presented 
earlier^® and compared with measurements for atoms 
in two-dimensional harmonic traps with a weak attractive 
interaction. In this work, it was shown that inhomogeneities 
must be considered for a correct determination of the quasi¬ 
particle effective mass. Here we discuss the role of inho¬ 
mogeneities in this experiment in more detail and propose 
ways to reduce their effect. We also show that, in the exper¬ 
iments of Ref. 24 made with ®Li atoms in three-dimensional 
harmonic traps, the inhomogeneity sets the line shape of 
the integrated rf intensity and should be considered for the 
precise experimental determination of the scattering length. 

The paper is organized as follows. In Sec. II A, we present 
the model and the calculation of ny(fc, t). The generic 
diagrams giving the momentum distribution are shown 
in Sec. IIB, and the leading contribution is evaluated in 
Sec. IIC. Without interaction in the final state, the analysis 
simplifies as shown in Sec. Ill, where we discuss the inter¬ 
play between the inhomogeneity and the Hartree shifts. We 
study final-state effects in Sec. IV: the effect of a finite life¬ 
time, the effect of inhomogeneous Hartree shifts, and other 
final-state effects which correspond diagrammatically to ver¬ 
tex corrections. Conclusions and perspectives are given in 
Sec. V. 


II. FINITE-TEMPERATURE, SECOND-ORDER RESPONSE 
THEORY FOR THE MOMENTUM DISTRIBUTION 

A. Description of the model 

The atoms are modeled as three-level systems with inter¬ 
nal states I a), a= 1,2,3. States |1) and |2) are assumed 
to interact most strongly, while state |3) has higher energy 
and will be the final state of the rf experiment; see Fig. 1. 
We consider that the three levels have the same dispersion; 
our results are readily generalized to the case where the 
dispersions are different in the initial and final states. These 
levels correspond, in practice, to atomic hyperfine states. 
The interaction between |1) and |2) is resonant and can be 
tuned by means of a Feshbach resonance.®® Once the field 
is set, the interactions between |1) and |3) and between 
12) and |3) are also set. Ideally, the latter interactions are 
small compared with the former. We consider hereafter 
fermionic atoms, and we assume translation invariance for 
simplicity. The extension to bosons is straightforward, and 
the formalism can be developed in real space if needed. Let 
be the creation operator for an atom in the state |a) 
with momentum Kk. The low-energy effective Hamiltonian 
is H = Hq + V , with 

Ho = Xi CD 

k 

^ ~ X ^2) 

kk'q 

al3 

We consider here the case of a local interaction between two 
atoms with center-of-mass momentum q. The detailed form 
of the interaction plays no role in our derivation, which is 



FIG. 1. (Color online) Rf spectroscopy of ultracold atoms. The 
hyperfine atomic levels |1), |2), and |3) are split by the Zeeman 
effect, such that a rf transition is allowed between states |2) and |3). 
The tunable many-body interaction between atoms of the cloud in 
states |1) and |2) shifts and broadens the level |2). The effects of 
this interaction are probed by comparing the frequency v of the rf 
transition with the frequency Vq of the noninteracting transition. 
Spin relaxation from |2) to |1) is forbidden by conservation of spin, 
meaning that the energy splitting between |1) and |2) is irrelevant. 
Levels 11) and |3) may be renormalized by interactions as well. 
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also valid for more general momentum-dependent interac¬ 
tions. For a contact interaction, the Pauli principle prevents 
atoms in the same internal state from interacting, and we 
can set = 0. 

The level separation Vq is typically in the 100 MHz range, 
and the rf radiation at this frequency has a wavelength of the 
order of meters. The rf pulse therefore induces momentum- 
conserving transitions. Let H'(t) be the time-dependent 
interaction between the rf radiation and the atoms, and let 
us assume that the allowed transition is between states |2) 
and |3). We have 

H'(t) = ^(t)^ + h.c.) =<f(t)2r^. (3) 

k k 

The function gives the time envelope and strength of 
the coupling. For later convenience, we define the operator 

rfc = rl = 4c2^ + h.c. (4) 


B. Generic diagrams for the momentum distribution 


We now expand the momentum distribution in the final 
state, Ufc = ^ 3 k^ 3 k’ powers of H'. In the grand-canonical 
ensemble, and in the interaction picture, we have 


{rikiO) =^rpnk[t), (5) 

with p = p _ ^ the chem¬ 

ical potential, and N the number operator. Since we work in 
equilibrium, the chemical potential p sets the populations 
of the three levels, and N is the total atom number. The 
evolution is given by n^Ct) = !7“^(t)nfc!7(t) with [/(t) the 
interaction part of the evolution operator. The zeroth-order 
term is obviously 

{nkit))^°^ = Trprik = {nk)H, ( 6 ) 


which gives the equilibrium thermal population of the final 
state. For a noninteracting system, this becomes 


(r^)c>H„ =/(efc+^'^'0-P). (7) 

where /(e) = -F 1)“^ is the Fermi distribution function. 
The usual setup is that the final state is initially empty, 
such that this contribution is negligible at low temperatures, 
fegT ^ hvg. The first-order term in H' is known, from 
standard linear-response theory, to be 


K(t))™ = . 


r f 


dt'{[nkitlHXt')])H 


r 00 

J —00 


dco 


k' 


At the second line, we have introduced (S’Cm), the Fourier 
transform of ifCt), and the equilibrium retarded correlation 


function of the operators uj. and pj./ in the system described 
by H. In the time domain, this correlation function is 


C,,(O = -^0(O(K(t),rv(o)])H, (8) 

with the time dependence of the operators governed by the 
evolution Because H, N, and Uk all conserve 

the number of atoms in the states |2) and |3), while the two 
terms in jk' do not, the correlation function (8) vanishes 
identically. 

The second-order response involves a double commutator 
and can be expressed in terms of a double-time retarded 
correlation function of the three operators uj., and Yk"' 




(2) = I 


r f 


r f 


de 


dt" 


J —00 J —00 


X {[[n,(t),H/tO],H'(t")])r 


dco dco' , ,,, 

- )tgl^co)Slico') 

2n 2n 

h'l-" 


(9) 


This contribution can be evaluated within the Matsubara for¬ 
malism. We find that the double-time correlation function in 
Eq. (9) is given by the anal 5 nic continuation to real frequen¬ 
cies of an imaginary-frequency function 
according to (see Appendix A) 


= 2n'Kk^T with integer n denote the even Matsubara 
frequencies. In the imaginary-time domain, the double-time 
function is defined as 

= (7’Tr^fc(T)rfc'(0)rfc"(T - tO>h, dD 

with the imaginary-time ordering operator, nj(.(T) = 
and similarly for Yk"('^ ~ '^0- The 
imaginary-time and imaginary-frequency functions are re¬ 
lated by 


rP 


^kYk'Yk' 


Xin,,in'j = 




( 12 ) 

The correlation function (11) is nonzero, because the two 
crossed terms in the product Jk'Tk" conserve the number of 
atoms of each flavor. Gathering these two terms, we get 


’^^krk'Tk"^'^’ ^ 

{T^clk[T)c^k^z)clkX0)C2kmclk,,[T - T0c3fc„(T - ^'))h 
+ {T,c]’^(t)C3^(t)c^,, ,(0)c3^,(0)c^'^„(t - T0c2fc„(T - tO)h- 

( 13 ) 
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FIG. 2. Diagrammatic representation of the double-imaginary 
time function (11). Diagrams (R) and (L) correspond to the first 
and second terms of Eq. (13), respectively, up to a minus sign 
(the correlation function equals minus the diagram, if standard 
diagrammatic rules are used). The two diagrams are identical, 
except for the direction of the arrows (R, right-handed, L, left- 
handed). The white dots denote the interaction of atoms with 
light, leading to transitions between internal states |2) and |3); the 
black dots denote the measured momentum distribution, and the 
hatched regions represent all interactions, including interactions 
with the state |1). 


The two terms can be represented by Feynman diagrams, as 
shown in Fig. 2. These diagrams have three entry points, 
one representing the measurement of the momentum dis¬ 
tribution and two representing the transitions induced by 
H' between states |2) and |3). Similar diagrams arise in the 
response theory of electron photoemission.This is to 
be contrasted with the bubble-type diagrams representing 
the transition rate Nf (see, e.g.. Ref 11). 


C. Leading contribution 



(I) (II) (lO 


FIG. 3. Leading term (I) and an example of vertex correction 
(II). In both cases, the two topologically inequivalent diagrams 
of kinds R and L must be considered. For a noninteracting final 
state, diagrams of type I' are the only nonvanishing contributions. 
Double lines denote the Green’s function, and single lines denote 
the noninteracting Green’s function in the final state. 


-(r.,c^^(T)c^^.(0))H. In imaginary frequency, they are 






1 

1 

- hvQ - SjCfc.imJ’ 


(14) 

(15) 


where ico^ = (2n-l- I'jnk^T, — p, and 

is the self-energy. Translating the two diagrams using the 

conventions of Fig. 2 gives 


+ - tO%(Ic, tO%(Ic, -t). (16) 

The minus sign associated with the fermion loop is canceled 
because the correlation function equals minus the diagram. 
We perform the Fourier transform in Eq. (12) using the 
spectral representation of the Green’s function, 

‘Sa{k,T)= d£Aa{k,e)-^- -, (17) 

J-OC 


We can distinguish two categories of diagrams, as illus¬ 
trated in Fig. 3. The justification for separating the diagrams 
of type I from “vertex corrections” of type II stems from the 
fact that, in usual experimental conditions, the interactions 
in the final state are small compared with the other interac¬ 
tions. If the former are exactly zero (V 13 = 1^23 = 1(33 = 0), 
all vertex corrections of type II disappear. Then, the two 
diagrams of type T (right- and left-handed) are the only 
nonvanishing terms, with the two propagators in state 3 
given by free propagators. Since diagram (T) can be derived 
from diagram (I) by taking the appropriate limit, we shall 
evaluate here diagram (I) and discuss the case of a nonin¬ 
teracting final state in the next section. Final-state effects 
are present in both type-I and type-II diagrams. We discuss 
final-state effects of type I in Secs. IV A and IV B and those 
of type II in Sec. IVC. 

The two diagrams of type I involve a single momentum, 
i.e., ^kk'^kk"- They can be expressed 

in terms of the Green’s functions for each atomic state. 
We define the fermionic Green’s function as (^„(fc,T) = 


where A„(fc,e) = —Im^„(fc,ia)n ^ e-l-i 0 ''')/ 7 r is the single¬ 
particle spectral function. This leads to 


Xy 

xiV — 

B ^ ico„ 


f 00 

dEde'de" 

J —00 

X A2{k, E^A^ik, E ^A^ik, e ) 


ico^ + — e' ico^ — in'^ — e" 


+ 


itUn — iHn — e' itOn 4- — e" 


(18) 


The frequency sums are evaluated in the usual manner^ ^ 
and yield terms proportional to either/(e), /(e 0 > or/(&")• 
For weak interactions in the final state, A 3 (fc, e) peaks near 
E = + Hvq. The terms proportional to / {s'}, and /(e") 

are therefore small at low temperature, like the zeroth-order 
term ( 6 ). We denote the contribution of the domi- 

nant terms proportional to /(e) and the contribution 
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of the other terms. We have 


nkYkYk^ riY 


dede'de" 


X A 2 {k, e)A 2 {k, e')Ao,{.k, e") 


a spectral function A 3 (fc, c) = 5(e — — hvo), the resolution 

function simplifies to 




1 


nt 2 

^i{e+hvyihgf^^,^ 

J —00 


(23) 


x/(e) 


ifln + e — e' — e + e" 
-1 1 


+ : 


iO.„ — e + e' iQ' + e — e" 

n _ 


(19) 


Making the analytic continuation as in Eq. (10), and insert¬ 
ing in Eq. (9), we obtain the leading contribution to the 
momentum distribution: 


(n,(t)) 0 ^) = 


dcA 2 (A:,e)/(e) 


/* 00 


dE'A^ik, e')^i^[e — e') 


( 20 ) 


The dimensionless function ^ accounts for the broadening 
effect due to the rf pulse: 


Equations (22) are one central result of this work. We 
use them to study the interplay of Hartree shifts and inho¬ 
mogeneity in two-dimensional (Sec. Ill B) and three- 
dimensional ^Li (Sec. Ill C) and to study the effect of a finite 
lifetime in the final state (Sec. IVA). 

The terms resulting from the frequency sum in Eq. (18), 
which have not been retained in Eq. (19), are 


YlkYkYk^ YlJ 


f 00 

dEdE'dE" 

^ —00 

y. A2{k,£')A2,^k,E'jA^ikjE ) 


fW) 


+ T 


.'O/ 


+ E — e' in'„ + E — e' I in^ + — e' + e" 


;cYi 


in„ — E + e' ifll — E + e' I + ifl'^ + e' — e" 


(24) 


J't(e) 


r“ dco 

I 271 hcj — e - I - iO "*" 

J —00 


— I 

Y 


r t 


dt'eY'-^^g[t'l 

) 

( 21 ) 


The main goal of rf spectroscopy is to determine the spectral 
function A 2 (fc,e). Eor weak interactions, this function is 
peaked near e = e^ — fj, = On the other hand, since 
the dispersions in the initial and final states are the same 
and only q = 0 transitions are possible, one expects to ob¬ 
serve, by varying the frequency v of the rf radiation, a signal 
peaking close to the frequency Vq of the noninteracting tran¬ 
sition. In order to make this more apparent, we introduce 
the detuning v = Vg ~ t', we change variables in Eq. (20), 
and rewrite it in the form 


{nfc(t)) 


(la) 


de^^'''^(e)A 2 (fc,?fc -l-hv - e) 

X fi^k+hv - e). ( 22 a) 


We have rearranged the terms by exchanging e' and e" in 
half of them. We proceed as above, and introduce again a 
resolution function: 


(rtfcCO)®'’’ = -/(?fc +(iVo) 


de<%™(e) 

X A 2 (fc, -b hv — e). (25a) 


We have pulled out a minus sign, because this term is nega¬ 
tive: It corresponds to a reduction of the thermal population 
in the final state as given by Eq. ( 6 ), induced by transi¬ 
tions to the initial state. These terms describe an inverse 
rf spectroscopy analogous to the inverse photoemission in 
condensed-matter systems. In the usual experimental prac¬ 
tice, they do not contribute because the atom cloud is pre¬ 
pared in a slightly out-of-equilibrium state, where level |3) 
is empty. The resolution function in Eq. (25) is 


This shows that the measured momentum distribution is the 
convolution of the occupied part of the spectral function 
with a dimensionless resolution function ^ . Under ideal 

conditions, the resolution function is proportional to 5(e), 
and the momentum distribution peaks near v = Vg, as ex¬ 
pected. The expression of the resolution function resulting 
from Eq. (20) is 




(* 00 


de'de' 


,fi^k+hv-E') 


/(?fc+K) 

X AgCfc, + hv - E')A2[k, + hv - e") 

r f ft' 


X Re 




dt e « 


SY) 


") 


(25b) 


Y\e) = 


dE'A^^k, e')^^[e + e' — 


. ( 22 b) 


It takes into account the renormalization of the final state 
by interactions, as well as the broadening due to the time 
envelope of the rf pulse. For a noninteracting final state with 


For a noninteracting final state, e' and e” are both equal 
to —hv. The upper limit of the t" integral can be extended 
from t' to t, correcting with a factor 1/2. This shows that 
Eq. (25b) reduces to Eq. (23) and that the two resolution 
functions are equal for a noninteracting final state. We 
finally note that, if t = -boo—i.e., if the measurement of the 
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momentum distribution is performed after the extinction of 
the rf pulse—the resolution functions are simply given by 

=—\S[E/h +27i:v)f, (26) 

ft 

with (S[cj) the Fourier transform of 


III. NONINTERACTING FINAL STATE 



e/{hAv) 


In this section, we neglect the interaction between the 
final state |3) and states |1) and |2) (V 13 = 1^23 = ^31 = 
V ^2 = 0 )- furthermore restrict to a short-range interac¬ 
tion such that = 0. The atoms are free fermions in the 
final state, and the nonzero matrix elements are V 12 = ^21 > 
describing the short-range interaction between the states 
|1) and |2). In this limit, diagram (F) in Fig. 3 gives the 
whole second-order response, and the momentum distribu¬ 
tion is the sum of Eqs. (7), (22), and (25). Because the 
spectral function in the final state is a 5 function, both res¬ 
olution functions are given by Eq. (23). In the context of 
electron photoemission, an analogous model known as the 
“sudden approximation” assumes a free-electron final state. 
In contrast to rf spectroscopy for cold atoms, however, this 
remains an approximation even in the ideal situation of a 
truly noninteracting final state, because other effects (in 
particular the surface) are usually neglected as well. 

We consider a monochromatic radiation of frequency v 
with a slowly varying envelope, such that the coupling in 
Eq. (3) is = <l(t)cos(27rvt). Assuming that the mo¬ 
mentum distribution is measured after the end of the pulse, 
and that the duration of the pulse is much longer than 1 /v, 
the resolution function is = |<#(e/ft)p/(4ft^) with 

<#(m) the Eourier transform of <?(t). For a square pulse of 
intensity Sq and duration Av“^ in the limit v » Av, the 
resolution function is 




7 7I<?o' 


’sin(7rx)' 

y ftAv ^ 

) 

nx 


e 

ftAv 


(27) 


For a Gaussian pulse of the same intensity at the maximum 
and a full width at half maximum Av“^ we have 


5?'''’(e) 



41n2 


exp 


(n^ 

In 4 


(28) 


The two functions are compared in Fig. 4. The validity of 
any approach based on equilibrium response is limited to a 
regime where the fraction of atoms transferred to the final 
state is small, or equivalently, the time t of the measure¬ 
ment is short compared with 2hjSQ, which is the period 
of Rabi oscillations between the states |2) and |3). Since, 
on the other hand, t > Av“^, this regime corresponds to 
<?o/(2ftAv) ^ 1. Our formalism is therefore valid as long 
as the amplitude of the resolution function is much smaller 
than unity. 


FIG. 4. (Color online) Resolution function for a rf pulse with a 
square (dashed line) and Gaussian (solid line) envelope in the 
case of a noninteracting final state. The width of the envelope in 
the time domain is fixed by Av“^ in both cases (full width at half 
maximum in the Gaussian case). The curves are normalized to the 
peak maximum for easier comparison. 


A. Free and nearly free fermions in a harmonic trap 

If all interactions are turned off, the spectral function is 
A 2 (k, e) = 5(e — For a homogeneous fermion gas, the 
momentum distribution is therefore simply 


{nfc) = <%"'’(ftv)/(^fc) 

-F [1-^'’'’(ftv)]/(^fc-Fftvo). (29) 


The first term corresponds to the atoms excited from the 
initial state, while the second term corresponds to the equi¬ 
librium thermal population of the final state, reduced by 
the transitions to the initial state. From here on, we assume 
that ftvg is large enough for the second term in Eq. (29) 
to be negligible. Equation (29) indicates that one can, in 
principle, determine the frequency Vq of the noninteracting 
transition, the resolution Av, and the temperature T by 
measuring the momentum distribution with all interactions 
suppressed: the energy-distribution curve (EDC) is just the 
resolution function, while the momentum-distribution curve 
(MDC) is controlled by the Fermi function. Experiments with 
homogeneous Fermi gases have not been conducted yet (for 
bosons, see Ref. 32). In this section, we study within EDA 
the modifications of Eq. (29) due to the nonhomogeneous 
distribution of atoms trapped in a harmonic potential, in 
two and three dimensions. The resulting equations provide 
a way of determining the total number of atoms, in addition 
to Vq, Av, and T, by fitting experimental EDCs and MDGs. 

Consider a harmonic trap described by the potential 
y(r) = (l/2)ma)V^. In dimension d, the number of atoms 
in state | 1 ) is related to the chemical potential by 


iVi = 




J 



d'^k 1 

(271)'' ^ ‘ 


(30) 


Eor free particles with a dispersion = ft^fc^/(2m), the 
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evaluation of the integrals gives 


B. EDC dispersion due to inhomogeneity and Hartree shifts 



where Li 2 and Lij are the di- and trilogarithm, respectively. 
Note that, unlike Eq. (31) suggests, p does depend on 
the particle mass m, because co^ = where k is 

the strength of the harmonic potential. To estimate the 
trap-averaged momentum distribution, we replace in 
Eq. (29) with / — p + V(r)), and we perform a spatial 

integration. The result is 


For free and nearly free fermions, the EDC line shape is 
not modified by the inhomogeneity, and the dispersion of 
the EDC peak tracks the intrinsic quasiparticle dispersion. 
However, if the self-energy depends on density, these con¬ 
venient properties are lost. In order to illustrate this in the 
simplest model, we consider the case of fermions subject to 
a short-range interaction, which is treated to lowest order, 
by keeping only the Hartree term. The momentum- and 
energy-independent Hartree self-energy in state |2) is given 
by 


/ 2Tck^T\^ f 

(r^fc)LDA = - I ^^2 J Li| l^-e ‘bt J (32) 

Note that (ujjjLDA is extensive and has the units of a nor¬ 
malization volume. Equation (32) shows that for free parti¬ 
cles in the EDA, the EDCs are not affected by the inhomo¬ 
geneities, because the latter do not change the energy of the 
|2) ^ |3) transition. The measured EDC line shape depends 
neither on the details of the density distribution in the trap 
nor on the momentum, but is entirely determined by the 
properties of the rf pulse. 

A qualitative understanding of the effects of interactions 
on the EDC and MDC curves may be gained by considering 
nearly free fermions. The simplest model is that of free 
fermions with an effective mass m*. With the caveat that 
such a model can only be envisioned as a low-energy ide¬ 
alization, this effective mass can be simulated by assuming 
for the bare fermions a self-energy. 



such that + = ft^fc^/(2m*) = E^., and the spectral 

function of the initial state is A 2 (k, s) = 5{e — — Sfc). Ne¬ 

glecting the population of the final state, the corresponding 
momentum distribution for a homogeneous gas is 

(nfc) = '%0')(hv-Sfc)/(Efc-/r). (34) 


S 2 = (g/iVo)ni, (35) 

with rti the density of atoms in state |1). The state |1) 
experiences a similar shift Ej = (g/iVo)n 2 . The dimension¬ 
less coupling g is positive (negative) for repulsive (attrac¬ 
tive) interaction, and Nq is the Fermi-level DOS, given by 
mfcp/(27r^ft^) and m/(27rft^) in three and two dimensions, 
respectively. The coupling g is related to the scattering 
length via g/iVo = ^nh^a^o/m and g = —l/ln(fcpa 2 D) in 3D 
and 2D, respectively. Equation (35) means that the energy 
of the transition is reduced (increased) with respect to the 
noninteracting value hvg for a repulsive (attractive) interac¬ 
tion. In a harmonic trap, the modification varies from the 
center to the periphery, and this contributes to a broadening 
and a momentum dependence of the EDC, resulting in a 
dispersion of the EDC peak, as we shall see. This disper¬ 
sion may by qualified “spurious”, because it is observed in a 
system where the transition does not actually disperse with 
momentum. 

With the Hartree term (35), the spectral function is 
A 2 (k, e) = 5(e — — E 2 ). For a homogeneous gas of den¬ 

sity Til, the chemical potential is set by the self-consistency 
condition 


d^k 

(^J_ 


deAi(fc,e)/(e) 


mkuT^ 2 


2 Tim 


M-(g/h^0)n2 


Lid —e ''b’' 


(36) 


The maximum of the EDC is at h(vo — v) = Sjt, and tracks 
the difference in the dispersions of the initial and final states. 
The dispersion of the EDC maximum is given by = 

(ft^fc^/2)(l/m* — 1/m). This means that, for m* > m, the 
peak moves towards lower values of the detuning Vq — v 
with increasing momentum k. 

Like for free fermions, the inhomogeneities due to trap¬ 
ping do not affect the EDC line shape for nearly free 
fermions, because the self-energy (33) does not depend 
on the local atom density. In such a gas, a plot of the 
quadratic EDC peak dispersion as a function of k gives the 
effective mass. The MDC profile also reflects the effective 
mass. Equation (31) gets corrected by a factor {m*/mY^^ 
because is defined in terms of the bare mass; in Eq. (32), 
the changes ^ and hv ^ hv — must be made in 
order to describe harmonically trapped nearly free fermions. 


and the momentum distribution (22) becomes 

K)=^0'^(hv-E2)/(?fc + E2). (37) 

In a harmonic trap, the local self-consistency condition reads 


ni(r) = - 


mk^T 

2ntY 


Lid —e 


(38) 


where p is fixed by the condition iV^ = J d'^rnyr').^^ The 
explicit expression for the momentum distribution in the 
harmonic trap is therefore 


(f^fc)LDA — 


dV ■ 


(hv -(g/iVo)ni(r)) 


exp I 






(39) 


+ l 
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FIG. 6. (Color online) Convergence of the measured EDC towards 
the intrinsic EDC (thin line) with improving the resolution by 
increasing the duration Av“^ of the Gaussian pulse. The EDCs are 
shown for g = —0.35 and k = 6 with the same parameters 

as in Fig. 5. 


FIG. 5. (Color online) Energy-distribution curves for harmoni¬ 
cally trapped atoms in two dimensions: effect of the Hartree 
shift. The EDC (39) for a Gaussian pulse with Av = 5 kHz are 
normalized and shown at different momenta (thick lines), for 
attractive (left) and repulsive (right) interaction. The thin lines 
show the intrinsic EDC corresponding to each curve [Av = 0, 
Eq. (40)], divided by two for clarity. The other parameters are 
co,/2n = 127 Hz, = 2000, and T = 100 nK. 


Figure 5 shows the EDCs calculated using Eq. (39), in a two- 
dimensional gas of atoms with Uj = n 2 = n, and using 
parameters typical for the experiment of Ref. 23. The max¬ 
imum of the EDC disperses towards lower (higher) values 
oivQ — v for attractive (repulsive) interaction. The width 
of the EDC varies with momentum and is larger than the 
expected resolution, which is Av\/21n4/7T = 3.12 kHz. 

The curves for attractive and repulsive interaction look 
similar in Eig. 5; however, the magnitudes of g are different. 
In fact, there is a systematic asymmetry between positive 
and negative g, because an attractive interaction tends to 
gather atoms near the center of the trap, leading to a more 
inhomogeneous density [see Fig. 7(a)]. The width of the 
EDC reflects the distribution of densities in the trap. This 
distribution is defined as D(n) = J d‘^r5(n — n[r)), and 
takes nonzero values for densities n between 0 and n(0). As 
shown in Appendix B, it is possible to rewrite the momentum 
distribution (39) as an integral over densities involving 
D(n) [Eq. (Bl)]. In dimension d = 2, the distribution D(n) 
can be evaluated explicitly (see Appendix B) . For an ideal 
resolution, the resulting momentum distribution is 


(r^fc>LDA « ^ 




0 ^ 

g JVo 


0 otherwise 


(40) 


with b(E) = — 1)“^. This is shown as the thin lines in 

Fig. 5 and corresponds to the Av ^ 0 limit of Eq. (39). The 
peculiar line shape (40), which depends on both the density 


distribution D(n) and the Fermi occupation factors, could 
be revealed experimentally by a moderate improvement of 
the resolution. Figure 6 shows the evolution of a typical line 
shape, as the full width at half maximum of the Gaussian rf 
pulse is increased from 0.2 to 1.0 ms. 

The dispersion of the EDC maximum with increasing 
momentum is plotted in Eig. 7(b) for various interaction 
strengths and temperatures. This “spurious” dispersion may 
be understood as follows. At each point in the trap, the 
minimum of the local energy band is the sum of the har¬ 
monic potential and the Hartree term. The k = 0 states 
are occupied throughout the trap and give contributions 
to the EDC with a Hartree shift ranging between 0 at the 
periphery and (g/iVo)n(0) at the center. The EDC for fc = 0 
extends therefore from v = 0 to v = (g/iVo)n(0)/h. As 
the momentum increases, the corresponding k states in the 
low-density regions at the periphery are above the chemi¬ 
cal potential, and their thermal population contributes less 
to the intrinsic EDC. The latter is depressed near v = 0 
and becomes asymmetric. Once filtered with the reso¬ 
lution function, the observed EDC disperses as seen in 
Eig. 5. We can be more quantitative in the limit T 0, 
where the intrinsic EDC (40) becomes a rectangular dis¬ 
tribution, constant for e^. ^ hv/g ^ n(0)/NQ, and zero 
otherwise. Convolved with the resolution function, this dis¬ 
tribution gives a peak whose maximum disperses quadrat- 
ically: E^QQ^^k) = (g/2)[e;j + n(0)/iVo]. This dispersion 
can be parametrized by a Hartree “effective mass” as 
^max(fc) - -EmaxCO) = (fi^fc^/2) (l/mn - 1/m). We then find 
m^ /m=l/(l-l-g/2). This is compared in Fig. 7 (c) with the 
mass calculated numerically for various temperatures. The 
density at the trap center can also be evaluated at T = 0: 
n(0)/iVo = HcJq y2iv570^+~^-^^ With this, we can calculate 
the full T = 0 dispersion, which is shown as thin lines in 
Fig. 7(b). 

At finite T, the dispersion is not quadratic, except close 
to fc = 0, and for g not too close to —1. As temperature 
increases, the particle cloud spreads more across the trap 
[see Fig. 7(a)], the density distribution D(n) sharpens, and 
the peak dispersion therefore gets weaker. The asymmetry 
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FIG. 7. (Color online) (a) Radial density, (b) dispersion of the EDC 
maximum as a function of momentum, (c) Hartree “effective mass” 
as a function of g, and (d) width of the EDC. In all graphs, dotted 
lines correspond to T = 150 nK, dashed lines to T = 100 nK, thick 
solid lines to T = 50 nK, and the thin solid lines give the analytical 
result for r = 0 (see text). In (d), the width is shown as a function 
of g for three momenta; the colors correspond to those used in 
Fig. 5. The model parameters are u>J2n = 127 Hz, iVj = 2000, 
and Av = 5 kHz. The minimum at g = 0 in (d) corresponds to the 
resolution of 3.12 kHz. 


between repulsive and attractive interaction is strongest at 
T = 0 and is reduced as temperature increases. The finite-T 
Hartree mass m^, deduced from the curvature of the dis¬ 
persion at fc = 0, is shown in Fig. 7(c). The relative mass is 
larger than unity for attractive interaction and smaller than 
unity for repulsive interaction. Its dependence on tempera¬ 
ture is linear for g > 0, but more complicated for g < 0; in 
particular, nonlinearities in the temperature dependence get 
stronger as g approaches —1, as can be seen in Figs. 7(b) 
and 7(c). 

As seen in Fig. 5, the EDC not only disperses due to inho¬ 
mogeneity, but also narrows with increasing momentum. At 
zero temperature, the width of the EDC has a complicated 
dependency on Av, which approaches a linear function 
of |g| as Av ^ 0, namely \glh{e^ — n(0)/iVo)|. At finite 
temperature and finite resolution, however, the EDC width 
behaves more like ~ g^, as shown in Fig. 7(d). At fc = 0, 
the width reflects the radial density distribution of Fig. 7(a): 
It is larger for attractive than for repulsive interaction of 
the same magnitudes and decreases with increasing tem¬ 
perature. The rounded behavior at g = 0 transforms into a 
linear behavior ~ \g\l\Jl + g as Av is reduced.At large 
momenta, the width is controlled by the Fermi edge rather 
than by D(n): It is resolution-limited at low temperature 
and increases with increasing T. 

The dispersion displayed in Figs. 5 and 7(b), as well as 



FIG. 8. (Color online) (a) Evolution of the Hartree “effective mass” 
with (u, for attractive (top curves) and repulsive (bottom curves) 
interaction and for three different temperatures, (b) Hartree “effec¬ 
tive mass” as a function of g for a r'^ trap and three temperatures. 
The thin solid line shows 1/(1 -I- g/2) for easier comparison with 
Fig. 7(c). (Inset) Radial density for two values of g and the same 
three temperatures. 


a width looking like ~ g^ as seen in Fig. 7(d), could easily 
be mistaken as a signature of dynamical interactions, since 
this is the expected behavior in a homogeneous Fermi liquid. 
Even the narrowing of the EDC with increasing fc could 
evoke the sharpening of quasiparticles when approaching 
the Eermi momentum. One feature, however, among these 
inhomogeneity-driven effects is contrary to the expected 
signature of interactions: The sharpening of the fc = 0 EDC 
with increasing T, due to the flattening of the atom cloud 
in the trap, cannot be reconciled with the expected increase 
of the scattering rate with T in an interacting system. 


An obvious way to reduce the spurious mass in ex¬ 
periments is to achieve a more homogeneous density. In 
three dimensions, this can be realized by means of a weaker 
confining harmonic potential. It does not work well in two 
dimensions, though, because the density distribution D(n) 
is flat at zero temperature: All densities between zero and 
the maximum density are equally represented, irrespective 
of the potential strength. This is illustrated by the fact that 
is given by 1/(1 -t- g/2) at zero temperature, which does 
not depend on < 1 )^. At finite T, D(n) does depend on n 
(see Appendix B), and thus m^/m approaches unity as 
is reduced, as shown in Eig. 8(a). Alternatively, one may 
obtain a more homogeneous density by means of a quartic 
trapping potential. Eigure 8(b) shows as a function of g 
when the harmonic trap is replaced by a trap. The latter 
was defined such that the potential equals the harmonic 
potential with = 127 Hz at a distance r = 10 pm. 

The inset shows that the density profile is flatter than in 
Eig. 7(a), and as a result is significantly reduced as 
compared to Eig. 7(c). The spurious dispersion can also 
be reduced, and in some cases even suppressed, by taking 
advantage of interactions in the final state; this is discussed 
in Sec. IVB. 
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C. Inhomogeneity and momentum-integrated rf intensity 


In this paragraph we briefly discuss the effect of Hartree 
shifts and inhomogeneity on the momentum-integrated rf 
intensity. In the particular case of a Gaussian density profile, 
which is a good approximation for three-dimensional gases 
at not too low temperature, and a Gaussian rf pulse, the 
trap-averaged integrated intensity takes a universal form 
depending on a single parameter oc a 3 nn( 0 )/(mAv). We 
compare this form with the measurements of Ref. 24, where 
Hartree shifts in ®Li mixtures were studied by rf spectroscopy 
as a means to determine the scattering length. 

The momentum integration of the rf intensity (39) yields 


■^LDA — 


(hv - (g/iVo)ni(r )). 


(41) 


For a balanced gas with = n 2 = n, this becomes a one¬ 
dimensional integral involving the density distribution; 

^ 00 

4da= dnD(n)n^^''^(hv-(g/iVo)n). (42) 

J —00 


We are interested in comparing this expression with the 
data of Ref. 24, which were obtained at T ~ O.yTp- At 
such temperatures, the density profile given by Eq. (38) 
is very close to a Gaussian in three dimensions: n(r) ^ 
n(0)exp{ — 71 [n(0)/iVi]^^^ r^}. For this Gaussian profile, 
the density distribution is 


D(n) = 


2JViy(ln[n(0)/n]/7r 

n(0jn 


0 ^ n ^ n(0) 


otherwise. 


(43) 


If the rf pulse has a Gaussian envelope, Eq. (42) becomes 


f grt(O) V 


4dA \ u a at ’ a 

y riAviVo Av 


f 2 fla 3 Dn( 0 ) V 


where the function is given by 


fOO 

71^ 9 

du Vue “ exp 
Jo 

- —(x-ae-“)2 

In 4 


(45) 


This function is displayed in Fig. 9(a). For a = 0, it is a 
Gaussian, and for a 7 ^ 0 it has an asymmetric shape with 
a cutoff at X a. Note that J^(a, x) = J^(—a, —x), so that 
symmetric curves are expected for attractive and repulsive 
interactions corresponding to identical values of the product 
WsdHO). 

Equation (44) can be fit to the ®Li data as shown in 
Fig. 9(b). The function captures remarkably well the 
peculiar asymmetry of the line shape in the presence of in¬ 
teractions (white dots), in contrast to the Gaussian used 
in Ref. 24, resulting in an excellent fit. We have set 
Av = 7.1 kHz, which corresponds to the Gaussian rf pulse of 
140 ps used in this experiment. For the data without inter¬ 
action (black dots), there is no further adjustable parameter. 




Vq — V (kHz) 

FIG. 9. (Color online) (a) The function J^(a, x) defined in Eq. (45) 
for various values of a. Note that the cutoff is at x ^ a. (b) Data 
of Ref. 24 (dots) and fits to Eq. (44) with Av fixed to 7.1 kHz 
(lines). A constant background was added to Eq. (44) for fitting. 


apart from the amplitude and a constant background. For 
the curve with interaction, the fit yields a = —7.15 ± 0.34. 
Considering that the average density is n ~ 2.4 x 10^^ cm“^, 
and that n = n(0)/2^^^ for a Gaussian profile, this corre¬ 
sponds to a scattering length Ujd = —35.4± 1.7 nm. The 
analysis of Ref. 24 yields a value consistent within the error 
bars, —31 ± 2.7 nm, but we believe that Eq. (44) provides a 
more accurate way of measuring the scattering length. 


IV DISCUSSION OF FINAL-STATE EFFECTS 


A. Resolution function for a Lorentzian final state 

We start this section by deriving the modifications due 
to the resolution functions (23) and (28), in the situation 
where interactions lead to a shift and a lifetime for the fi¬ 
nal state of the rf transition. We introduce these effects by 
means of a phenomenological self-energy S 3 — in the 
final state, where S 3 is the energy shift and Fj is the scatter¬ 
ing rate. These two quantities are, in principle, related by 
causality and should be of the same order of magnitude for 
weak interactions. The corresponding lifetime of the final 













11 


state is T 3 = and the spectral function reads 


A3,{k,e) 


_£ 3/71 _ 

i£-^k-hvo-T.3f + rl' 


(46) 


The noninteracting result (23) gets modified like this: 


rt 2 

dt' 

— 00 

(47) 

The overall magnitude of the resolution function vanishes on 
time scales larger than T 3 , because atoms in the final state 
decay. Besides, the energy dependence of the resolution 
function is also affected. In order to find out how, we 
perform the time integration explicitly for the case of a 
Gaussian pulse of full width at half maximum Av“^. In the 
relevant limit Av ^ v, the formula replacing Eq. (28) is 


>(Ia) 


(e) = 


A1 2 1 / fl 




X exp 


71 

exp 

r 

16 In 2 

ln4 V 

■ 


In 4 

y hAv J 


2r3t 


1 + erf 2 


!Vhi2A 


£ + So — ir - 


vt ■ 


2Vin2 hAv 


(48) 



(f- +S3 )/(?iAv) 



t (Av“^) t (Av“^) 


This complicated expression has an interesting time de¬ 
pendence (Fig. 10). The resolution function is even and 
centered at the energy s = —S3, and it is significantly 
non-Gaussian when the time delay t of the measurement— 
counted in Eq. (48) from the maximum of the pulse 
envelope—is comparable to the width of the pulse. For large 
times t » Av“^, the erf function approaches one, and the 
energy dependence of Eq. (48) measured from e = —S 3 
is identical to the noninteracting result (28), shown in 
Fig. 10(a) as “Fourier limited”. The width of the resolu¬ 
tion function takes off for measurement times of the order 
of Av“^ and increases roughly linearly with decreasing t 
[Fig. 10(b)]. The peak intensity of is largest shortly 

after the pulse maximum and decreases for longer times 
[Fig. 10(c)]. 

Measurements can be done in the regime where the res¬ 
olution function is Fourier limited, provided that the time 
t is smaller than the lifetime T 3 , but sufficiently large, that 
the real part of the argument in the erf function is large and 
positive. These requirements read 


1 

161n2 T 3 AV 


Av"^ 


< t < T 3 . 


(49) 


Clearly, such a regime does not exist unless T 3 » 
Av“^/(4\/ln2) or r 3 ^ (\/ln 2 / 7 r)hAv, as illustrated in 
Fig. 10(c). 


B. Hartree shifts in the final state 

In Sec. Ill B, we assumed = n 2 , which is justified if the 
fraction of atoms transferred to the final state is small. A 


FIG. 10. (Color online) (a) Energy dependence of the resolution 
function (48) for a final-state scattering rate 03 = hAv, correspond¬ 
ing to a lifetime T 3 = Av“^/( 47 i), and for increasing measurement 
times (broader to narrower). The time t is measured from the 
maximum of the Gaussian pulse envelope, such that t = Av“^ 
corresponds to a measurement time one full width after the pulse 
maximum, (b) Full width at half maximum of the resolution 
function relative to the Fourier limited value and (c) maximum 
intensity as a function of the measurement time and scattering 
rate. The arrows pointing to the left in (c) indicate the time T 3 , 
and those pointing to the right the time ^ Av~^. 


more accurate modeling of experiments on balanced gases 
would be to take n 2 = Uj — n^, where is the number of 
atoms in the final state. If is a fraction / of tij, we may 
write n 2 = (1 —/juj and = frii. Let us furthermore take 
into account the interactions and g 23 between states 
|1) and |3) and states |2) and |3), respectively, in addition 
to the interaction gi 2 (which was denoted g in Sec. IIIB). 
Treating all interactions at first order, we find that the level 
I a), a = 1,2,3, is shifted by the self-energy 

Sa = (g:/iVo)ni, (50) 

with g] = gi 2 (l - /) + gis/, g 2 = gi 2 + g 23 f, and g^ = 
gi 3 + g 23 (l ~ The shift S 3 of the final state is larger at 
the center of the trap than at the periphery and will there¬ 
fore contribute to the spurious mass m^. We assume that 
I S 3 1 remains much smaller than hv^, such that interaction 
effects related to the thermal population of the final state 
are negligible. 

The resolution function reflects the shift of the final state: 
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+ S 3 ). As a result, Eq. (39) is replaced 

with 


(f^Jc)LDA — 


d'^r- 


exp 


(hv - ^^ni(r)j 

r ft-M+lmM^r^+Cg^/NojniCr) 

I, 


(51) 


+ 1 


One sees that the width of the rf signal is now controlled by 
g* —g* instead of gi 2 - If the parameters (interactions and/or 
transferred fraction /) can be arranged such that gl = g^, 
then the dispersion of the final state locally follows the 
dispersion of the initial states, and no spurious dispersion 
should be observed. 

An explicit expression for the Hartree “effective mass” 
in the presence of final-state shifts can be derived in two 
dimensions: The ideal momentum-distribution line shape 
(40) is replaced with 



lo 


otherwise. 


(52) 

In the limit T ^ 0, this becomes again a steplike distribu¬ 
tion, whose center disperses quadratically with momentum. 
Proceeding as in Sec. Ill B, we find 


l+gt-g; 

m l + gt-ig*2 + g;)/2’ 


(53) 


which is indeed unity if g 2 = gl- 


C. Vertex corrections at low temperature and density 

The vertex corrections of type II describe final-state effects 
going beyond the self-energy renormalizations of the final 
state. We estimate such effects in this section and indicate 
how they could be implemented to improve the theoretical 
description of rf measurements. In the context of electron 
photoemission, specific vertex corrections were shown to 
describe the production of plasmons^® or phonons^^ during 
the photoexcitation process. These phenomena are not rel¬ 
evant for cold-atom systems, but other interesting effects 
take place, related to the spatial correlations among the 
dilute atoms. We proceed in two steps, in order to iden¬ 
tify the important vertex diagrams. First, we consider the 
regime k^T hvg and eliminate all diagrams that require a 
thermal population of the final state. Then we organize the 
remaining diagrams according to the number of hole lines 
in the initial states in the spirit of the low-density expansion 
for the self-energy.^'* We furthermore assume a short-range 
potential, such that the interactions are blocked by the 
Pauli principle. 

This analysis, outlined in Appendix C, shows that the 
most important vertex diagrams are those represented in 
Fig. 11. Diagram (II.Rl) describes the correlated state of 




FIG. 11. Dominant vertex corrections in the low-temperature 
low-density limit. The shaded boxes represent a particle-particle 
ladder series (pseudopotential). Ip and Ih stand for a particle or a 
hole in state |1), respectively, and similarly for the other states. The 
diagrams give a significant contribution only if the time ordering 
of the various vertices is such that all lines marked as particles 
go to the right (see Appendix C). The diagrams shown are right 
handed; there are two equivalent left-handed terms. 


three atoms during the rf conversion. Before the conversion, 
the atom in state | 2 ) is entangled with an atom in state | 1 ). 
This entanglement is preserved once the atom |2) is excited 
by the rf radiation to state |3). If the interaction V 32 
attractive, this process enhances the effect of the final-state 
interaction V 33 by increasing the probability that the excited 
atom has an atom in state |1) nearby. If V 32 repulsive, this 
process keeps the excited atom away from atoms in state 
|1), reducing the effect of V 13 . The effect of the final-state 
interaction V 23 > on the other hand, is limited by the Pauli 
principle: Just after the conversion, the atom |3) sits in the 
correlation hole of the former atom | 2 ) and is kept away 
from other atoms in state |2). The exclusion principle indeed 
forbids any contribution like (II.Rl), where the atom |1) 
would be replaced by an atom |2). The converted atom can 
nevertheless interact with atoms in state | 2 ), either directly 
(self-energy corrections of the 3p line) or via the exchange 
process represented by the diagram (II.R2). In this process, 
the converted atom interacts with an atom in state | 2 ) that 
is present above the Fermi energy, such that the interaction 
does not produce a new hole. The atom |2) eventually 
recombines with the hole left by the conversion, while the 
atom 13) is converted back to an atom |2) above the Fermi 
sea. 

In the experimental setup of Ref. 23, the scattering length 
ai 2 measuring the interaction V 12 between states | 1 ) (|F = 
9/2, mp = -9/2)) and |2) (|F = 9/2, = -7/2)) is close 

to a Feshbach resonance and was tuned on the attractive 
side from 0 to —300, in units of the Bohr radius Ug. The 
interactions V 13 and V 23 between states |1) and |3) (|F = 
9/2, mp = —5/2)) and |2) and |3) are both nonresonant 
and repulsive and correspond to scattering lengths 013 = 
-t-250aB and 033 = -l-130aB. In this configuration, we expect 
that the effect of ^13 is enhanced by the attractive V 12 in 
the correction (II.Rl), while V 23 onfy contributes through 
the exchange process (II.R2). We speculate that most of 
the extra broadening observed in the measurements, with 
respect to the theory including V 12 , but neglecting final-state 
interactions,^^ is the result of these processes. A self-energy 
broadening due to the direct interaction between |2) and |3) 
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in the initial state (accounted for in the type-I diagram of 
Fig. 3) is unlikely because such contributions require at least 
two holes in the final state and are suppressed by a factor 
Q-hvo/kgT Appendix C). If the number of atoms excited 
in the final state is not too small, the self-energy in the final 
state associated with both and V 23 rnay also induce, in 
addition to energy shifts at lowest order, some broadening 
of type I, which enters the resolution function, as shown in 
Sec. IV A. Explicit evaluations of the vertex corrections in 
Fig. 11 and of other self-energy effects are left for future 
works. It will be interesting to see whether and how these 
final-state effects change the line shape of the rf signal. 


Appendix A: Anal5n:ic continuation of second-order response 
functions 

In this appendix, we show that the second-order retarded 
susceptibility, defined in terms of the double commutator 
in Eq. (9), corresponds by analytical continuation to the 
imaginary-time correlator (11). We switch to a slightly 
lighter notation, set ft = 1 , and compute the second-order 
change of the expectation value of an observable A, in the 
presence of a perturbation H' = BF(t), where B is an ob¬ 
servable and F(t) is a classical field. The second-order 
correction is 


V CONCLUSION 

We have presented a theoretical description of the rf spec¬ 
troscopy of cold-atom systems, based on the second-order 
response theory at finite temperature. The difference be¬ 
tween the usual golden-rule approach and this new descrip¬ 
tion is that the latter focuses on the number Nf of atoms 
transferred to the final state, while the former focuses on 
the transition rate Nf. The second-order response approach 
accounts for the finite energy resolution implied by the enve¬ 
lope and the finite duration of the rf pulse and allows one to 
classify the various contributions using Feynman diagrams. 
The issue of inhomogeneity represents a challenge for the 
interpretation of rf experiments performed on interacting 
Fermi systems. Due to the density dependence of the self¬ 
energy, the rf line shape varies across the cloud. We have 
studied this effect at leading order in the density within the 
EDA and found that the static local Hartree shifts induce an 
apparent dispersion of the rf signal, similar to the disper¬ 
sion expected in a homogeneous interacting Fermi gas from 
dynamical effects of higher order in the density. For three- 
dimensional gases with a Gaussian density profile, we have 
derived a simple expression for the momentum-integrated 
rf line shape, which takes into account the finite resolution 
and the inhomogeneous Hartree shifts. 

Final-state effects are another challenge for rf experi¬ 
ments. We have considered the simplest of them, resulting 
either from a lifetime or from the interplay of inhomogeneity 
and Hartree shifts in the final state. More subtle final-state 
effects, such as those resulting from the spatial correlations 
between atoms, are described by vertex corrections. We 
have proposed a scheme to classify these terms and identi¬ 
fied those which dominate at low temperature and low den¬ 
sity. A numerical evaluation of the corresponding diagrams 
is needed to tell whether these effects change significantly 
the line shape of the rf signal. 


r t 


{A(t))® = (-1)2 


r 


dti 


dtn 


J —00 J —00 


X {[[A(t),H'(ti)],H'(t 2 )])„. (Al) 


The ensemble average is taken over the eigenstates of the 
time-independent Hamiltonian H, such that invariance by 
translation in time applies: {[[A(t),H'(ti)],H'(t 2 )])H = 
{[[A(t - ti),H'( 0 )],H'(t 2 - Using this, and the 

expression of H', we can write 






■ti,t-t 2 )F(ti)F(t 2 ), (A2) 


with the second-order susceptibility defined as 

X ([[A(t),B(0)],B(t-tO])H. (A3) 

Introducing the Fourier transform of the various quantities 
in the integrand of Eq. (A2) leads to the analogous of the 
second line in Eq. (9) : 


(ACt))® = 


f 00 T T / 

dco dco 


2n 2n 




^ Z®('^>'^OB(m)F(mO. (A4) 


Out task is to show that 

^ co + iO+,in' ^ co' + iO+), (A5) 


where iO') is the Eourier transform of the 

imaginary-time correlator 
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Eor this purpose, we show that the spectral representations 
of the functions and are identi¬ 

cal. 

(2^ 

Let us start with ■ Splitting the imaginary-time inte- 
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grals to take into account the time ordering, we have 




rP 

0 

rP 


dTdT'e'^^e'^'^V®(T,TO 


dTe"' 


dr' e'^ (A(t)B(t — tOB(O))^ 


rP 


+ 


dr' e'^ ^ {A(r)B[0)B[r - tO)h 


(A7) 


To perform the time integrations, we introduce a complete 
set of eigenstates of H, H\a) = BqIu), we use the expression 
of the thermal average, ((■ • ■)) = Z“^Tre“^^(- ■ •we insert 
two times the identity U = |a) (a|, and we use the expres¬ 
sion of the imaginary-time operators, e.g., A{r) = . 

The averages in the square brackets of (A7) become 


(A(t)B(t - tOB(0)>h 

^ abc 

(A(t)B(0)B(t - tO>h 

= -2(a|A|b)(b|B|c){c|B|a>e^“^^'-«e-^»^e^'(^-^'^. 


The r and r' integrations in (A7) are now elementary and 
yield, after making use of the property e'^P = e'^P = 1, 


‘^j^\in,in')=jY,(a\A\b)(b\B\c}{c\B\a) 

abc 


iQ -|- -|- Efj 


Q-PEa _ g-PSc 


iQ -f E„ 


+ 


e-PEt _ q-Pe, 

in — Eu+Er 


Q-PEa _ g^-PEc 


+ 777 


g-PEb _ ^-PE, 


in' + E„ 


+ 777 


in — Eh+Er 


■ (A8) 


A similar calculation leads to the spectral representation of 
the real-time susceptibility. We start from 


The four terms of the double commutator are expressed as 

{A(t)B(0)B(t-tO)H 

= i2(a|A|b)(b|B|c){c|B|a)e-^^“e'^“''e-'^'>'e'^'^'-''^ 
(B(t-tOA(t)B(0))„ 

= -2(a|A|b)(b|B|c){c|B|a)e-^Ae;£„t'e-<£btgi£,(t-t') 

^ abc 

(B(0)A(t)B(t-tO)„ 

= i2(a|A|b)(b|B|c){c|B|a)e-^^'e'^Ae-‘'^..t'e'7(t'-t) 

^ Tbt 

(B(t-tOB(0)A(t))H 

= i2(a|A|b)(b|B|c){c|B|a)e-^^'>e'^“'e-'^'>''e'Ah'-t). 


We perform the time integrations in (A9) with the help of 
the identity 


p 00 

dte'“'B(t) = 

J—00 

and obtain, using the notations co 
(X) H“ iO^, 


i 

CO + I0+’ 

■^ = CO + lO"^ and co 


/+ 


%^\co,co')= -^(a|A|b)(b|B|c){c|B|a) 

abc 

I f g-PEa _ g-PEc g-pEt _ ^-PE, \ 

co^ + co'^ + Bjj — Bj, y co'^ + E^ — E^ co'^ — Ef, + y 

By exchanging the dummy variables tj and t 2 in the expres¬ 
sion (A2), we see that the susceptibility (A3) can also be 
defined with the arguments t and t' exchanged. We could 
therefore use an alternate definition of the susceptibility, 
which shows explicitly the symmetry under the exchange of 
the time arguments, e.g., t') + t)] instead 

of Eq. (A3). Exchanging the time arguments in Eq. (A3) is 
equivalent to exchanging the two frequencies co and co' in 
Eq. (A9). After performing this symmetrization, we obtain 
the alternate definition of the susceptibility: 


X^\co,co') = ■^i^(a|A|b)(b|B|c){c|B|a) 

^ abc 

I f g-PEa _ g-pE, g-PE, _ ^-PE, 

^ m++ + B„ -Bf, l^m+TB^-B, m+-Bb + B, 

g-PEa _ g-PEc Q-PEt _ g-PE, \ 

H-TT-1- -r - ) . (AlO) 

CO + Bq — B^ CO — Efj E^ J 

Comparison of Eqs. (A8) and (AlO) proves (A5). 


(2)/" /\ 
X)jjXco,co) 


p 00 

dtdt'e'“'e'“'''0(t)0(t'-t) 

J —00 

x([[7l(t),B(0)],B(t-O])H. (A9) 


Appendix B: Momentum density and density distribution 

By inverting the analog of Eq. (38) for n 2 (r), one ob¬ 
tains an expression for p — (l/2)mcu^r^ — (g/iVo)ni(r) as 
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a function of n 2 (r). Inserting this expression into Eq. (39) 
gives 


{rik) 


LDA • 


d^r ■ 


-(g/iVo)ni(r)) 




_ f 2rf h 
V niksT J 


Nd/2 


n2ir) 


where Li“^ is the inverse of the polylogarithm function. If 
ni(r) = n 2 (r) = n(r), the r-dependence of the integrand 
stems from n(r), and the spatial integration can be con¬ 
verted into a density integration, by introducing the density 
distribution D(n) = J d'^r 5 [n — n[r)): 


('^fcZtDA — 


^ D(n)i% 0 'Hfiv-(g/iVo)n) 
dn-;=- , ■ (Bl) 


1 -e 


I \inksTj 


For an ideal resolution, 5(^s), we have simply 

DihvNo/g) 


(r^)c>LDA « 




_ f 2rf '\'^I2 ^yjV^ 
ymkgT J g 


This expression can be made more explicit in dimension 
d = 2. On the one hand, Lq ^(x) = 1 — e“^, and on the other 
hand, the density distribution can be evaluated explicitly. 
We have 


D(n) = ^ 


2nro 

\n'[ro)\ 


0 


0 ^ n ^ n(0) 
otherwise, 


where n'[r) is the derivative of the radial density n(r), n(0) 
is the density at the trap center, and n(ro) = n. Differentiat¬ 
ing Eq. (38) with respect to r, one finds 


2Tir 271 

|n'(r)| tuco^Nq 


1^- 2^oj^r^-(g/No)n(r) 

1 + g + e ‘bJ- 


The exponential in the square brackets can be expressed as 
a function of n(r) only, by inverting Eq. (38) as above. For 
r = to, on thus gets 


27rro 271 r f ^ ^ 

|n'(ro)| mm^iVo 


where b(s) = — 1). The resulting expression for 

the momentum distribution in two dimensions, and for an 
ideal resolution, is given in Eq. (40). Interestingly, the 
functional dependence of the density distribution on n, and 
consequently the dependence of the momentum distribution 
(40) on V, does not involve the total particle number Mj; 
only the cutoff depends on iV^ via n(0). 


Appendix C: Classification of vertex corrections 

The upper line in the diagram of type R in Fig. 2 corre¬ 
sponds to a hole in the final state |3), as implied by the 



FIG. 12. Contributions to the right-handed (R) and left-handed 
(L) diagrams of zeroth order in the interaction. The vertices are 
ordered horizontally by increasing imaginary time from left to right. 
2p and 3p indicate particle lines in states |2) and |3), respectively, 
while 2h and 3h indicate hole lines. 


ordering of the times, e.g., r — t' < r. The lower line 
corresponds to a particle in the final state. Conversely, in 
the diagram of type L, the lower line corresponds to a hole 
(0 < t) and the upper line to a particle. In both cases, the 
vertical line describes either a particle or a hole in the initial 
state 12), depending upon the ordering of the times r — t' 
and 0. This is illustrated in Fig. 12 in the case of type-I 
diagrams. Each hole in the state |3) entails an occupation 
factor /(fis) ~ which is negligible if the thermal 

population of the final state is negligible. One such factor 
is canceled if—and only if—the time r' — r can reach the 
value (This applies to R diagrams; the same statement 
with r' — r replaced with r — r' applies to L diagrams.) The 
reason is as follows. The Green’s function for a free hole 
propagating between times Tj and T 2 is 
whereas for a free particle it is All time 

dependencies from the various particle and hole lines in 
state 13) cancel, except at the two conversion vertices (o), 
leaving only the dependence \ Upon performing 

the time integrations as specified by Eq. (12), a factor 
is generated if the time r' — r (t — t' for L diagrams) is 
allowed to reach the value 13. This explains the behaviors 
indicated in Fig. 12. In all cases, there is one hole in the 
final state (3h line)—hence a factor —that is canceled 

for right-handed diagrams if t — t' < 0 and for left-handed 
ones if T — t' > 0. The two types of contributions were 
denoted (la) and (Ib) in Sec. IIC. 

Since the cancellation of the final-state hole occupation 
factor can only work once, we conclude that any diagram 
with more than one hole in the state |3) carries at least one 
factor and is exponentially small if k^T Hvq. In 

particular, all corrections of the density vertex (•) imply a 
connection between the lines 3p and 3h that cuts the 3h line 
and thus contains at least two holes in the final state. The 
first-order corrections of the conversion vertices (o) which 
survive in the limit fcgT ^ hvg are displayed in Fig. 13. 

At higher orders in the interaction, we classify the vertex 
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corrections like in the low-density expansion of the self- 
energy.^"^ A self-energy diagram containing p hole lines, for 
instance a particle-hole ladder at order p -I-1, is propor¬ 
tional to Since p ^ — oo as the density n—* 0 at any 
finite temperature, the contributions with one single hole 
dominate in this limit. These contributions are given by the 
particle-particle ladder series. Similarly, the vertex correc¬ 
tions with one single hole in either of the initial states |1) 
or 12) are expected to dominate at low density. Figure 11 
shows the two contributions which we consider as the most 
important vertex corrections at low density. Both contain 
a single hole in the final state, and a single hole in one of 
the initial states. Any further decoration of these diagrams 
with interaction lines introduces new hole lines. The two 
right-handed first-order terms of Fig. 13 may be obtained 
from the diagram (II.R2) by removing one of the interaction 
boxes and evaluating at first order. 


(R) 


(L) 






FIG. 13. First-order right-handed (R) and left-handed (L) con¬ 
tributions to the o-vertex corrections, which survive in the limit 
fegT C hvg. The imaginary times t — t', 0, t", and t are ordered 
horizontally as in Fig. 12. Any modification in the ordering of 
times produces at least one factor . 
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